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Abstract 

vo 

£\| , The paper focuses on unconditionally optimal error analysis of the fully discrete Galerkin 

finite element methods for a general nonlinear parabolic system in R d with d = 2,3. In 
terms of a corresponding time-discrete system of PDEs as proposed in [22], we split the 
error function into two parts, one from the temporal discretization and one the spatial 
discretization. We prove that the latter is r-independent and the numerical solution is 



< 

* 1 n 

bounded in the L°° and W '°° norms by the inverse inequalities. With the boundedness of 



the numerical solution, optimal error estimates can be obtained unconditionally in a routine 
way. Several numerical examples in two and three dimensional spaces are given to support 
our theoretical analysis. 

Key words: Optimal error estimates, unconditional stability, Galerkin, nonlinear parabolic 
system 



^ ■ 1 Introduction 

CO 

There are several numerical approximations schemes in the time direction for the numerical 
solution of nonlinear parabolic equations (systems). Linearized (semi)-implicit schemes are the 
most popular ones since, at each time step, the schemes only require the solution of a linear 
system. However, time-step size restriction condition is always a key issue in analysis and 
computation. For many nonlinear parabolic systems, error analysis of finite element methods 
(or finite difference method) with linearized semi-implicit schemes in the time direction often 
requires certain time-step conditions. See [TJ [THl [201 [23 [29] for the Navier-Stokes equations, 
[HI [40] for the nonlinear Joule heating problems, [HI [THl 021 032 GS] for flows in porous media, 
[El EE] for viscoelastic fluid flow, [2S1 [39] for the KdV equations, [27] for the Ginzburg- 
Landau equations, [3l [33j H2] for the nonlinear Schrodinger equations and [THl [EH E7J for 
some other equations. Such time-step size restrictions may result in the use of an unnecessarily 
small time-step size and extremely time-consuming in practical computations. To study the 
error estimate of linearized (semi)-implicit schemes, the boundedness of numerical solution (or 
error function) in the L°° norm or a stronger norm is often required. If a priori estimate for the 
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numerical solution in such a norm cannot be provided, one may employ the induction method 
with an inverse inequality to bound the numerical solution, such as 

\\U% - R h u n \\ L oo < Ch- d / 2 \\U% - R h u n \\ L 2 < ChT d l 2 {r m + h r+1 ), (1.1) 

where U£ is the finite element solution, u is the exact solution and Rh is certain projection 
operator. A time-step size restriction arises immediately from the above inequality, particularly 
for problems in three dimensional spaces. Most previous works follow this idea. A new approach 
for unconditionally optimal error analysis of a linearized Galerkin FEM was pesented in our 
recent work [22] , also see [24], where the error function is split into two parts, the spatially 
discrete error and the temporally discrete error, 

IK - a? || < IK - u n \\ + \\u n - u%\\, (i.2) 

where U n is the solution of a corresponding time-discrete parabolic equations (or elliptic equa- 
tions). Optimal estimates for the second term can be obtained unconditionally in a traditional 
way if suitable regularity of the solution of the time-discrete system can be proved. More re- 
cently, unconditionally optimal error estimates were established for a nonlinear equation from 
incompressible miscible flow in porous media. In [221 123] . analysis waw given only for a linear 
FEM and a low-order Galerkin-mixed FEM, respectively. 

In this paper, we consider a general nonlinear parabolic equation (or system) 

— -V-(a(u)Vu)=g(u,Vu,x,t) (1.3) 

in a bounded and smooth domain in M. d (d = 2 or 3) with the boundary condition 

u = on dn (1.4) 

and the initial condition 

u(x, 0) = uq(x) for i£!l, (1-5) 

where g 6 C 2 (IR) is a general nonlinear source. The general equation is of stronger nonlinearity 
than those in [22\ [23] and many physical equations are included. We apply linearized backward 
Euler Galerkin method with r-order finite element approximation (r > 1) for the general non- 
linear system. We focus our attention on the unconditional convergence (stability) and optimal 
error estimates of the linearized Galerkin FEMs. A key to our analysis is the a priori estimate 
of the numerical solution. We apply the splitting technique proposed in [22j [23] to bound the 
numerical solution UJf in L°°-norm and W 1,00 -norm, such as 

< \\RhU n \\ + \\UZ-R h U n \\ 

< \\R h U n \\ w x,oo + Ch- d ' 2 \\UZ - R h U n \\ H i 

<C + Ch- d / 2 h k (1.6) 

where k > d/2. Then with the boundedness of ||?7?||^i,oo, optimal error estimates can be easily 
established unconditionally in the routine way of FEM error analysis. 

The paper is organized as follows. In Section 2, we present linearized backward Euler 
Galerkin FEMs for the general nonlinear parabolic equations f|1.3[) - (|1.5j) and introduce our no- 
tations. In Section 3, we prove the boundedness of the numerical solution in the W 1,co norm 
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in terms of a corresponding time-discrete system, in which a rigorous analysis on the regularity 
of the solution to the time-discrete PDEs is given. Due to the boundedness of the numerical 
solution in the W 1 ' 00 norm, we present unconditionally optimal error estimates in Section 4 in a 
simple and routine way. Numerical examples in two and three-dimensional spaces are presented 
in Section 5. Numerical results confirm our theoretical analysis and show that no time-step 
condition is needed. 

2 Fully discrete Galerkin FEMs 

Let 7r/i be a regular division of 17 into triangles Tj, j = 1, • • • ,M, in M 2 or tetrahedras in M 3 , 
and let h = maxi<j<jv/{diamTj} denote the mesh size. For a triangle Tj with two corners (or a 
tetrahedra with three corners) on the boundary, we let Tj denote the triangle with one curved 
side (or a tetrahedra with one curved face). For an interior triangle, we simply set Tj as Tj itself. 
Finite element spaces on {Tj} have been well defined, e.g., see |32l I41j . For a given triangular 
(or tetrahedral) division of 17, we define the finite element space 

Vh = {vh G C(flh) '■ Vh\Tj is a polynomial of degree r and v h = on dfth} 

so that Vh is a subspace of H^Tlh). Let G : 17^ — > 17 be a coordinate transformation such that 
both G and G _1 are Lipschitz continuous and, for each triangle Tj, G maps Tj one-to-one onto 
Tj gl]. We define an operator Q : L 2 {VL h ) -> L 2 (17) by Qv(x) = v(G~ 1 (x)) for x G 17. Then we 
set 

V h = {Gv h : v h G Vh}. 
Easy to see that Vh is a finite element subspace of Hq(Q) and 

1 1 to — il/jto||LP < C||to|| H /r+i, P (Q)/i r+1 , for r > 1, 2 < p < oo, 

where Tlh = G^hG 1 and Tlh : Co(17ft) — >• Vh is the Lagrangian interpolation operator of degree 
r. 

Let = io<£i <••• < t]v = T be a uniform partition of the time interval [0,T] with 
t n = tlt and let u n = u(x,t n ) for n = 0, 1, • • • , N. For a sequence of functions {/ n }^L , we 
define 

D T f n+1 = J - J —, for n = 0,l,--- ,JV-1. (2.1) 

r 

A simple linearized backward Euler Galerkin method for the problem (|1.3j> — (|1.5j) is to seek 
U™ +1 eV h ,0<n<N-l, such that 

[D T U^\ v) + {a(U r h l )VU^\Vv) = {g(U r h l , W£, x, t), v) (2.2) 

for any v G Vh, with the initial condition U® = Hh u o f° r r > 2 and U® = R\uq for r = 1, where 
R\ is a projection operator defined in Section 3.2. 

With a linear approximation to the nonlinear source term, an alternative linearized scheme 
is defined by 

[p T Vj+\v) + (a(Uj:)VUj: +1 ,Vv) = (gS,v) + ( 9 ^ +1 ,v) + (g n 2 ■ VTJf\v), V v G V h , 
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where g% = g(U%yU%,x,t), g? = V l9 (U%, VU%, x, t) and g$ = V 2 g(U^VU^,x,t), with V l9 
and V25 denoting the gradient of g with respect to the components U and VfJ, respectively. 
The corresponding linearized Crank-Nicolson schemes can be defined similarly with classical 
extrapolations [12]. 

In this paper, we only focus our attention on the linearized scheme (|2.2p . The analysis 
presented in this paper can be extended to the second linearized scheme and many other schemes. 
We assume that g GC 2 (R x W 1, x £2 x [0, T]) and a £ C 2 (M) satisfies the weak ellipticity condition 

a(s) > for s£R. (2.3) 

3 Boundedness of the numerical solution 

In this section, we assume that the solution to the problem (|1.3p - (|1.4p exists and satisfies that 

IMU°°((0,T);ff3) + l|dHU°°((0,T);£ 2 ) + ll<9HU 2 ((0,T);ff 2 ) + \\<htu\\l?((p,T);L*) < M ; (3- 1 ) 
for some positive constant M, and we prove the following theorem. 

Theorem 3.1 Suppose that the system U.3\) - [T75\) has a unique solution u satisfying the reg- 
ularity condition §3. 1}) . Then there exist positive constants C and ho, independent of n and h, 
such that the finite element system A2.fy) admits a unique solution {Uf^}^ =1 when h < ho, and 

\\UhU°° + ll VC/ h < C - ( 3 - 2 ) 

To prove Theorem 13. 1( we introduce a corresponding time-discrete equation as proposed in 
[221 [23]: 

D T U n+1 - V • {a{U n )VU n+1 ) = g(U n ,VU n ,x,t n ), (3.3) 

with the boundary condition U n+l = on dQ. and the initial condition U = uo- 

In the following two subsections, we estimate the error functions u n — U n and U n — UJ^, 



respectively, where U n is the solution of the time-discrete system (13.31) . 

For the simplicity of notations, we denote by C a generic positive constant and by e a generic 
small positive constant, which depend solely upon M, £1, T, a and g, and independent of r, h 
and n. 

3.1 The time-discrete solution 

By the regularity assumption ()3.ip . we have u G W l, °° . We set 

K = ||n||ioc.(nx[o,T]) + \W u \\L°°(nx[Q,T}) + 2 • 
Q K = [-K,K] d+1 x nx [0,T]. 

Then, by the regularity assumptions on g and a and the ellipticity condition (|2.3p . there exist 
positive constants ok and Ck such that for \s\ < K and (a,/3,x,t) G Qk, 

o~k < o-(s) < C K , 

\a'(s)\ + \a"(s)\<C K , (3.4) 

d d 

\g(a,P,x,t)\ + \d a g{a,j3,x,t)\ + ^ \dp.g(a,P,x,t)\ + ^ \d Xj g(a, f3, x, t)\ < C K . 

j=i 3=1 
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Lemma 3.1 (i^'-estimate of elliptic equations [9]) Suppose that v is a solution of the boundary 
value problem 

Av = /, in 0,, 
v = 0, on dfl, 

where Q E M. d , d = 2,3, is a smooth and bounded domain. Then 

\\v\\ H i<C\\f\\ H i-2, 1 = 2,3. (3.5) 

In this subsection, we explore the regularity of the solution to the time-discrete system (|3,3|) 
and present an error estimate for u n — U n . 

Theorem 3.2 Suppose that the system \1.3\ )- [T75\) has a unique solution u satisfying ( pOP - 
Then the time-discrete system \3. 3]) admits a unique solution {U n }^ =0 such that 

N 

max \\U n \\ H3 + V T\\D T U n \\ 2 H 2 < C , (3.6) 

n=l 

and 



max \\u n - U n \\ t7i < Cqt (-i>7\ 
o<n<7v" l|jH ~ ^-'J 



where Cq is a positive constant independent of n, h and r. 



Proof For the given U n , (13. 3p can be viewed as a linear elliptic boundary value problem. With 
the first inequality in (|3.4p and classical theory of elliptic PDEs, the equation (|3.3p admits a 
unique solution U n+l in H . Let e n = u n — U n . Here we only prove the estimates fj3.6[) - (|3. 7|> . 
First, we prove by mathematical induction the inequality 

\\U n \\ L ~> + ||W n || L sc < for n = 0,1, ••• ,iV (3.8) 

under the condition r < r for some positive constant r . Since U° = uq, the inequality (|3.8p 
holds for n = 0. Now we assume that the inequality holds for < n < k. 

Let e n = u n — U n . From (|1.3p - ()1.5p and (j3.3j) . we see that e n+1 satisfies the equation 

D T e n+1 - V • (a(U n )Ve n+1 ) (3.9) 
= R n+l + V • ([fj(u n ) - a(U n )]Vu n+1 ) + 5 (n n , Vn n , x, t) - g{U n , VU n ,x, t), 

with the boundary condition e n+1 = on dQ and the initial condition e° = 0, where 

R n + l = d t u n+1 - D T u n+1 + V • {{a{u n ) - a{u n+1 ))Vu n+1 ] 
+g(u n ,Vu n ,x, t) - g(u n+1 ,Vu n+ \x, t) 

is the truncation error due to the time discretization. By the regularity assumption (|3.1[) . we 
have 

TV 

max ||i2 n || L 2 < C, V T\\R n \\ 2 r2 < Ct 2 . (3.10) 

~ ~ n=l 



5 



Multiplying the equation (|3.9p by e n+1 , we obtain 

Dr(l\\e n+1 \\l?) +^\\Ve n+1 \\ 2 L2 < ({a{U n )-<?{u n ))Vu n+1 , Ve" +1 ) 

+ ( g(u n , Vu n , x, t) - g(U n ,VU n , x, t),e n+1 ) +(R n+ \ e n+1 ) 



By (13. 41) . we have further 

\g(u n ,Vu n ,x,t) - g(U n ,VU n ,x,t)\ < C K (\e n \ + |Ve"|), 
\a(U n )-a(u n )\<C K \e n \. 

It follows that 

^(^ll^ +1 |li 2 )+^l|Ve" +1 ||i 2 

< e||Ve n+1 ||2 2 + Ce^\\e n+1 \\l 2 + e\\Ve n \\ 2 L2 + Ce- l \\e n \\\ 2 + Ce^WR^ 1 ]^. 
By choosing e < gk /4 and using Gronwall's inequality, there exists t\ > such that when r < n 



i 



||e n+1 || L 2 <C^r||i? m+1 ||| 2 y<CT, 0<n<fc, (3.11) 

^m=0 ' 

which implies that 

\\U n+1 \\ L2 < \\u n+l \\ L2 + ||e" +1 || i2 < C, (3.12) 
\\D T U n+1 \\ L 2 < \\D T u n+1 \\ L 2 + \\D T e n+1 \\ L 2 <C. (3.13) 

Applying Lemma 13.11 for the linear elliptic equation (13. 3b with the induction assumption gives 
the H 2 estimate 

\\U n+1 \\ H2 < C\\D T U n+1 \\ L2 +C\\Va(U n ) - VU n+1 \\ L2 +C\\g(U n ,VU n ,x,t)\\ L2 

< C||VC/ n || L oo||V[/ n+1 || L2 + C < C, (3.14) 

for < n < k. By the Sobolev interpolation inequality, 

||e fc+1 ||L~ < C\\e k+l \\ l - 2 d/ % k+l \\%t < Cr l - d l\ (3.15) 
Again we multiply the equation f)3.9|) by — Ae n+1 to get 

D T Q|| Ve n+1 ||| 2 ) + (a(U n )Ae n+1 , Ae n+1 ) 

< Ct~ l \\V<j{U n ) ■ Ve n+1 ||| 2 + Ce- 1 (|| J R n+1 ||| 2 + ||V • {{a{u n ) - a{U n )\Vu n+l )\\ 2 L2 
+ \\g(u n ,Vu n , x, t) - g(U n ,VU n , x, t)\\ 2 L2 ) + e|| Ae" +1 1|| 2 . (3.16) 

By (|3.4|) . the Sobolev interpolation inequality and the induction assumption, we have 

\\Va(U n ) ■ Ve n+1 \\ L2 < C||VC/ n || L oo||Ve n+1 || L2 < C||Ve n+1 || L2 , 



(u n ,Vu n ,x,t) -g(U n , VU n ,x,t)\\v < C\\e r 



\H 1 
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and 

\\V-[(a(u n )-a(U n ))Vu n+1 ]\\ L , 

< C||e n+1 || H1 ||V^ +1 || LO o +C||e n+1 || L3 ||An™ +1 || Le 

<C||e" +1 || Hl . 

Using Lemma |3. II and choosing a small e, the inequality ()3.16j) reduces to 

^r(||Ve n+1 ||| 2 ^ + ||e n+1 ||^ 2 < C\\e n+1 \\ 2 H1 + C\\e n \\ 2 H1 +C\\R n+1 \\ 2 L 2 . 
With GronwalPs inequality, we see that there exists a positive constant t 2 such that when r < t 2 , 



I n+1,,2 

\ e Wh 1 



+ X]r||e m+1 ||^ 2 < Ct 2 (3.17) 



m=0 

which together with (|3.ip leads to 

||e" +1 |Ui <C 2 r, ||C/" +1 || H2 <C 2 , 
\\D T U n+1 \\ m < C 2 , El=o4 D rU m+1 \\ 2 H 2 

for < n < k. 



(3.18) 



Moreover, we rewrite the equation (|3.3p as 

- AU n+1 = — !— ( 5 (C/ n , W n , x, t n ) - D T U n+1 + X7a(U n ) • W n+1 ) . (3.19) 
a{U n ) 

By Lemma O and (ET3D . 

||^ +1 || H 3 < ||— !— {g(U n ,VU n ,x,t n ) - D T U n+1 + Va(U n ) -VU n+1 ) ||^ 

< C\\U n \\ m +C\\g(U n ,VU n ,x,t n )\\ m + \\D T U n+1 \\ m +C\\Va(U n )-VU n+1 \\ H1 

< C + CII^H^ + C||C/ n || if 2||VC/ n+1 |U 00 +C||V^ n || L oo||C/™ +1 || H2 

< C + C\\VU n+1 \\ L °o 

<C + e||?7 n+1 || H3 + Ce^WU^ 1 ^, (3.20) 
which in turn implies that 

\\U n+1 \\ Wi <C 3 (3.21) 
if we choose e < 1/2. By the Sobolev interpolation inequality, 

||Ve fc+1 || LOO < C||e fc+1 ||^ /4 ||e fc+1 ||^ 4 < Cr 1 "^ 4 

which with fj3. 15f) shows that there exists T3 > such that 

\\U k+l \\ L ~ + ||VC/ fc+1 || Loo < ||n fc+1 || Loo + \\Vu k+1 \\ L ov + ||e fc+1 || LO o + || Ve k+1 \\ L °o < K (3.22) 

for r < T3. Thus (|3.8p holds for n = k + 1 when t < tq := minjri, t 2 , T3} and the induction is 
closed. Prom ([3Tllfj) - (|33TT ). we see that (j3lU)-(|3l7]) hold when r < r . 
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Secondly, we prove that (|3.6[) - (|3.7|) hold for r > To- We assume that maxi< n <fc ||[/ n ||#3 < 7^ 
for some positive constant 7^ (which may depend upon to) since ||[/°||#3 = H^ollija < C. From 
()3.3p . it is easy to see that 

k I 

\\U k+1 \\v < ( £ C\\g(U n , VU n ,x, t n )\\\A 2 <C lk . 

Then we apply Lemma 13.11 to (|3.19p . Via a similar approach as (|3.20p , we can derive that 

\V nJrX \ H z < C lk := 7 fe+ i. 
Since N = T/t < T/tq, we take C4 = max <A;<Ar 7fc so that 

\\U n+1 \\ H3 < C 4 , 0<n<N-l 

which further shows that 

\\D T U n+1 \\ HS <C a t~\ 

\\u n+1 - U n+1 \\ H i < ||n n+1 ||^i + \\U n+1 \\ H i <M + C 4 . 

Thus the induction is complete for t > tq. 

Combining the two cases, we complete the proof of Theorem 13.21 I 

3.2 The proof of Theorem [HQ 

For n > 0, let R^ +1 : H^(Q) ->• V h be a projection defined by 

(a(U n )V{w - Rl +1 w),Vv) = (3.23) 

for all v G Vh, and we set := R h . With the regularity of U n proved in Theorem 13.21 we have 
the following inequalities: 

jjn+l _ R n+l un +l^ L2 + qun+1 _ 1|^ < C\\U n+1 \\ H 3h 3 , if T > 2 (3.24) 

jjn+l _ R n+l un +l^ L6 + h \pn+l _ Rn+ljjn+1 1| ^ < C\\U n+1 \\ W 2, 6 h 2 , if r = 1 (3.25) 

K +1 U n+1 \\ w x^ <C||[/" +1 || w i,oc, (3.26) 

D T R r h l+1 U n+1 \\ w i,e < C\\D T U n+1 \\ w t,e, (3.27) 

D T {U n+1 - K n h +l U n+1 )\\ H ^ < Ch 3 , for r > 2, (3.28) 

D T (U n+1 - Rl +l U n+1 )\\ L 2 < Ch 2 , for r = 1, (3.29) 

where (|3.24[) - (j3.25p are standard error estimates of elliptic equations, (|3.26p - (|3.27p follow from 
|28j and the references therein, (|3.28[) - (j3.29[) can be proved in a similar way as in [221 123] . 
The following inverse inequalities will also be used in our proof. 

d_d 

IMIlp < Chp i \\v v G Vh, 1 < q < V < oo, 
||Vu||lp < C/i _1 ||v||lp, v G Vh, 1 < p < oo. 

Let 

Ki= max ||C/ n || w i,oo + max \\R h l U n \\ w \, ao + 2, r* = min{r, 2} . 

0<n<N 0<n<N 
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By the regularity assumptions for a and g, there exist o* K and C* Kx > such that 

\a(s)\ + \a'(s)\ + \a"(s)\ + \g{a, f3,x,t)\ + \d a g(a, f3,x,t)\ 
+ Y,j=i\ d Pj9(<x,P,x,t)\ + Ej=i I dad/^a, 

+ T, d j=i\d Xj dp.g{a^,x,t)\ + Ei,j=i\d^g(a, f],x,t)\ < C* Kl , (3.30) 
<r(s)>a* Kl , (3.31) 

for all s £ [— Ki,K{\ and (a,(3,x,t) 6 Qxi- 

Proof of Theorem \3.1\ We shall prove 

Ml~ + ||Vt^|| L oo < K x , n = 0,1, (3.32) 

n 

ll^|li 2 + E r ll Ve ™Hi^ 5 o /i2r * +2 ' ( 3 - 33 ) 

m=0 

simultaneously by mathematical induction, where = R^U n — UV'. It is easy to see that 
the inequalities (|3.32p - (|3.33p hold for n = 0. So we can assume that (|3.32j) - ()3.33|) hold for 
< n < k. By (|3.30p - (|3,3ip . the coefficient matrix of the linear system (|2.2p is symmetric 
and positive definite. Therefore, the system (j2.2|) admits a unique solution U£ +1 in Vh for 
< n < N- 1. 

Since the solution of the time-discrete equation (13. 3p C/ n satisfies 



(D T C/" +1 , W ) + {a(U n )VR n h +1 U n+1 ,Vv) = (g(U n ,VU n ,x,t),v), Vu 6 V^, 

it follows that e^ +1 satisfies the equation 

{D T el + \v) + (a(Uj;)Ve n h + \Vv) (3.34) 
= -((o-(U n ) - a(UJi))VRl +1 U n+l ,Vv) + (D T (R r h l+l U n+l - U n+l ),v) 
+ (g(U n ,VU n ,x,t) -g(UZ,VUZ,x,t),v), V « e V h . 

Now we estimate the last three terms in the above equation, respectively. For the first two 
terms, we see that 

| {(a(U n ) - a(UJi))VR n h +1 U n+ \ Vv) | < C(\\e n h \\ L 2 + h r ' +1 )\\Vv\\ L 2 

and 

\(D T {Rl +l U n+1 - U n+1 ),v) \ < \\D T (R^ +1 U n+1 - U n+1 )\\ H ~i \\v\\ H i . 
We rewrite the third term by 

{g(U n , VU n , x, t) - g(UJi,VUf:, x, t),v) (3.35) 
= (g(U n ,VU%,x,t)-g(Uj:,VUj:,x,t),v) + (g(U n , VU n ,x,t) - g(U n , VU%,x,t),v) . 

Here, we have 

\[g(U n ,VUZ,x,t)-g(UZ,VUZ,x,t),v)\ < C(|K|| L2 + //* +1 )|M| L2 . 
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By Taylor's formula, we get 

g(U n ,VU n ,x,t)-g(U n ,VUZ,x,t) 

= V(U n -UJi)- ( V 2 g(U n ,(l - s)VU n + sVU%,x,t)ds 
Jo 

= V(U n -£/£)• ( [V 2 g(U n , (1 - s)VU n + sVU%, x, t) - V 2 g(U n ,VU n , x, t)} ds 
Jo 

+ V • ({U n - UZ)V 2 g(U n ,VU n ,x,t)) - (U n - C#)V • (V 2 g(U n , W", x, t)) 
where V 2 g denotes the gradient of g with respect to the second conponent. Therefore, 



g(U n , VU n , x, t) - g(U n ,VUJt, x, t), *j 

< \\v{u n - uz)\\l 12/6 \\v\\ L6 + - uzilvWvWv 



+ 



(U n -UJi)V 2 g(U n ,VU n ,x,t),Vv 



< C(\\VeUU /5 + h r * +1 )\\v\\ Hl + C(\\e%\\ L2 + h r ^)\\v\\ m 



Substituting v = e£ +1 into (flTMD . we obtain 



n+l||2 



L 2 ) + CTA'-JlVe 



n+1 1|2 
h Hi 2 



(3.36) 



^ ^llY7^ n +l II 2 i 11^ _i /-< ^ — J- II „n. II i n ^— 1 1| w^, 

^ e ll Ve /i II x-2 + eHe^ || L 2+o 6 e ||e h || L 2 + 676 || Ve 



,n+l i|2 



lll„nl 



,-l|IV7 "||4 

'ML 12 / 5 



+ Ce- 1 ||D T (^ +1 C/ 



.n+ljjn+l _ jjn+l\\\2 



.i +Ce- l h 



li,2r*+2 



By an inverse inequality and the induction assumption, 

[|Ve£[|£ u/B < CCo^-^llVe^lH, < ^||Ve^||| 2 
if CC h 3 - d/3 < 1. With this estimate, we take h < ecr^/^CV) and sum up (|3.36p to get 



(3.37) 



n-l 



\ p n+l\\2 

\ e h IIl 2 



m=0 



By (explicit) Gronwall's inequality, we derive that 

n 

||„n+l||2 , ^||V7 m + 1 ll 2 <r r h 2r * 



2r*+2 



*+2 



(3.38) 



m=0 



for < n < k. 

To complete the mathematical induction, we need to prove (|3.32|) - (|3.33|) for n = k + 1. For 
this purpose, we consider two cases. 

Case I: r > 2. In this case, r* = 2 and we can apply inverse inequalities for (|3.38p to get 

||Ve£ +1 || LOO < Ch- d l 2 - l \\e n h +1 \\ L 2 < C 9 h r *- 3 / 2 , 
\\el+%~<Ch- d ^e? +1 \\r,*<C»h r *- l /\ 



\ Hi 2 
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which implies that 

\\U% +1 \\ w i^ < \\e n h +1 \\ wltOB + \\Rl +l U^ +l \\ w ^ < Kt (3.39) 

when Cgh 1 / 2 < 1. This completes the induction for r > 2, and (|3.32p - (|3.33p hold for all 
0<n<iV-landO<r<T. 

Case II: r = 1. In this case, r* = 1. To get the boundedness of ||Ve^ +1 ||L°°, we present the 
-H^-estimate with an additional induction assumption: 

\\Vel\\ L , + \\D T el\\ L ,<C l h 2 . (3.40) 

From the initial condition, we see that (|3.40p holds for n = and we can assume that it holds 
for < n < k. We substitute v = -D r e^ +1 into (|3.34p . With a similar approach to (|3.35p . we 
obtain 

^(iv^JVe^l^+p^lli, (3.41) 

< [D T a(UZWet, VeJI) - ((<r(CP) - <r(^))Vi^ +1 £/"+\ £ T Ve£ +1 ) 
+ C(p T (^ +1 ^ +1 - C/" +1 )|| L2 + ||Ve£||| 4 + \\e n h \\ L , + h 2 )\\D T e n h +1 \\ L2 
+ ((E7* - U%)V 2 g(U n ,VU n , x, t) , D T Ve 



n+l\ 
h ) ■ 

Using (|3.26p and the Sobolev embedding inequalities 

\\D T R%U n \\ L °° < C\\D T RZU n \\ w x, e , \\D T U n \\ w i,e < C\\D T U n \\ H 2. 
The first two terms of the right-hand side of the equation (|3.4ip are bounded by 
(D T a(K)VelVe^ < \\D T UZ\\ L <~ || Ve£||| 2 

<(p T e^||^ + ||D r it!^||^)||Ve%||| 2 
< C(h- 3 / 2 \\D T e n h \\ L i + \\D T U n \\ H ,)\\Ve n h \\ 2 L2 

and 

- {(a(U n ) -a{UZ))VRl +1 U n+1 ,D T Ve n h +1 ) 
= -D T ({a{U n ) -a(UZ))VRl +1 U n+ \Ve n h +1 ) 
+ (D T [(a{U n ) - a(UZ))VRl+ l U^\Vel) 

< -D T ((a(U n ) - a(UJi))VRl +1 U n+ \Ve n h +l ) 

+ C(\\D T (U n - UJt)\\ L2 + \\(U n - UZ)D T U^ L2 )\\VRW n h~\\Ve% +1 \\ L2 

+ c\\u n -un L 4DrVK +1 u n+1 ]\\ L 4^e n h +1 \\ L2 

< -D T ((a(U n ) - a(U%))VRl +l U n+ \Ve n h +1 ) 

+ e\\D T e n h \\h + Ce-\\\D T U n+1 \\h + 1) II Ve^ 1 1|| 2 + e\\e n h f H1 + Ch\ 
Moreover, we have 

[(U n - Ufi)V 2 g(U n ,VU n ,x, t) , D T Ve n h +l ) 
= D T {(U n - U n )V 2 g(U n ,VU n ,x,t),Ve n h +1 ) 
- {D T [(U n - UZ)V 2 g(U n ,VU n ,x,t)],Ve n h ) 
< D T {(U n - UZ)V 2 g(U n ,VU n ,x,t),Ve n h +1 ) + e(\\D T e n h \\ 2 L2 + \\D T (U n - R n h U n )\\ 2 L2 ) 
+ Ce~\l + \\D T U n \\ H i)\\Ve n h \\h + <H\\h + ll^ n " Eft||£e) • 
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With the above estimates, the inequality (|3.4ip reduces to 



D T 



n+l 



a(C/")Ve 



L 2 



+ \\D T e 



n+l i|2 



< Ce-\l + p r ^||^)(||V^||l 2 + ||Ve™ +1 ||| 2 ) + e(\\D T e n h \\h + KWh) 
+ Ce- 1 \\D T (iq+ 1 U r,+1 - U n+1 )\\ 2 L2 + Ce^WD^R^U 71 - U n )\\ 2 L2 
+ D T J n+1 + Ch\ 



where 



J n + l = - ((a(U n ) - <j{U%))VRl +l U n+ \ Ve™ +1 ) 
+ {(U n - U n )V 2 g(U n ,VU n , x, t) , Ve™ +1 ) . 

Furthermore, summing up (|3.42p gives 



n 

< Cn T i 1 + \\ D rU m+1 \\ 2 H 2) ||Ve™ +1 ||i 2 + Ch 4 



m=0 



where we have noted that 



Since 



|J n+1 | < e\\Ve n h +1 f L2 + Ce- l h\ 



Y,r\\D T U m+1 \\ 2 H2 <C, 

m=0 

by Gronwall's inequality and (|3.6p . (|3.43p further reduces to 

n 

||Y7n+li|2 , 11 n _m+li|2 / (~i h 4 



m=0 



for < n < k, provided r < T5 for some positive constant T5. 
Now by an inverse inequality, we have the estimate 



\e h Ww 1 - 00 — ll e h II// 1 ^ ^13" , 



and so 



7-rn+lii ll^'^+J-ll 1 II U 71 ^ 1 TT n +± II ^ Z<^ 

Illy 1 ' 00 — ll e h IIVK 1 ' 00 + 1 1 -"-ft U HVF 1 ' 00 — A l 



„n+l 



n+lj-rn+li 



(3.42) 



(3.43) 



when C13/1 1 / 2 < 1. It suffices to choose Cq > 1+Cs and C\ > l+2\/Ci2 so that the mathematical 
induction is closed when t < T5. It follows that (|3.32p . (|3.33p and (|3.40p hold for all n = 1, 2, iV 
when r < T5. 

When r > rs and r = 1, we can see from (|3.38p that 



U2 



(3.44) 
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which together with an inverse inequality implies that 

\K +1 \\ W ^ < Ch- d / 2 \\e n h +1 \\ m < C^h 1 ' 2 < 1 
ifh< C? b /t|. Therefore, 

+ ||V*7£ +1 || L ~ 

< \\R n h +1 U n+1 \\ L °o + \\VR n h +1 U n+1 \\ L oo + ||e^ +1 || L - + ||Ve" +1 || L oo < K\. 

It suffices to choose Co > 1 + C$ and C\ > 1 + Cut^ 1 so that the mathematical induction 
is closed for r > r 5 . It follows that (j!T32]) . (l3T33]l and (^IQ|) hold for all n = 1,2, ...,iV when 

Combining the two cases, we complete the proof of Theorem 13,11 I 

Remark 3.1 We have proved Theorem [37T] for any r-order Galerkin FEMs under the regularity 
assumption ()3. 1 1) . Based on the classical theory of finite element approximation and interpo- 
lation, this assumption is enough to obtain optimal error estimates for linear and quadratic 
Galerkin FEMs. In fact, the optimal I? error bounds for r = 1 and r = 2 have been given in 
(|3.33p . Since the estimates in (|3.33p are r-independent, by an inverse inequality, 

\\e n+l \\m < Ch r . 

By Theorem 13.21 and the projection error estimates in (|3.24p . we have optimal error estimates 
for the linear and quadratic Galerkin FEMs, which are summarized below. 

Corollary 3.1 Under the assumptions of Theorem \3.1\ there exist positive constants C and 
ho such that when h < ho, 

\\Uh ~u n \\ L 2 < C(T + h r+1 ) (3.45) 
\\U% -u n \\ m <C(T + h r ) (3.46) 

for r = 1 or r = 2. 



4 Error analysis 

Based on the boundedness of the numerical solution proved in the last section, one can easily 
obtain optimal error estimates of any r-order Galerkin FEMs under corresponding regularity 
assumptions, by following the classical approach of FEM analysis. Also it is possible to present 
the optimal error estimate for as we did in Section 3.2. However, this requires a rigorous 
analysis for stronger regularity of the time-discrete system. For simplicity, we follow the classical 
FEM approach and give a brief proof of optimal error estimates of the fully discrete Galerkin 
FEM. In this section, we assume that the solution to the initial-boundary value problem (11.31) - 
(jl.4p exists, satisfying (|3.ip and the following condition 

\\ U \\L^((0,T);H r + 1 ) + l|dH| L 2(( 0iT ) ;// max(r,2)) < C. (4.1) 

Let 9% = U% — RfrU 71 where R 7 ^ is the elliptic projection defined by 

(a(u n )V{w - Rl +l w), Vv) = (4.2) 

for all v <E Vh, and we set R° h := r\. Easy to see that (|3.24p - (|3.29p also hold for the projection 
operator R^ + . 
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Theorem 4.1 Suppose that the problem \ 1.3\) - [T75\) has a unique solution u satisfying h3. 1\) 
and Ii4-1\ ) f or some positive integer r. Then there exists he > such that for h < h§, 

\M -u n \\ L 2 <C(r + h r+1 ). (4.3) 

Proof Note that the error function 8 r £ satisfies the following equation: 

{D T 9l+\v) + [a{UZ)Vei + \Vv) (4.4) 
= -{(a(u n ) - a(UH))VR n h +l u n+ \Vv) + {D T {R n h +l u n+1 - u n+1 ),v) 
+ (g(u n ,Vu n ,x,t) - g(U%,VU%,x,t),v) + (R n+1 ,v) for all v E V h , 

where R n+1 is the truncation error satisfying (|3.1U[) , 

By the same approach as used in the proof of Theorem 13. 1\ we can derive that 

\((a(u n )-a(Uj:))VR n h +1 u n+1 ,Vv)\ < C\\0Z\\„\\Vv\\„, 

\{D T (R n h +1 u n+l -u n+l ),v)\ < \\DAR n h +1 u n+1 -u n+1 )\\ L2 \\v\\ L2 , 
| (g(u n , Vu n ,x, t) - g(U%,VU%, x, t),v)\ 

< C||W£|| l2 ||W£|MMIl 2 + Ch r+1 \\v\\ m + C(||^|| i2 + h r+1 )\\v\\ m . 
Since ||V#^||lo° < C as implied by Theorem 13.11 by taking v = in ([4.4p . we derive that 

^ii«0+^ii v oii 2 

< e(liv^ +1 ||i 2 + ||v^ +1 || i2 ) + crHWWh + \K +1 \\h) 

+ Ce- 1 \\D T (R n h +1 u n+1 - u n+1 )\\ 2 H ^ + C(\\R n+1 \\h + h * r+2 ) ■ 
By Gronwall's inequality, there exists a positive constant tq such that when r < tq, we have 

N 

max \\6%\\h + y r||V^||| 2 < C 17 {r 2 + h 2r+2 ). 

0<n<N n L ^ ll A 111, '\ ) 

Therefore, (|4.3p holds when r < tq. 

For r > r 6 , by TheoremEU \\U% - u n \\ L 2 < C < Ct^{t + h r+1 ). 
The proof of Theorem 14.11 is complete. I 
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5 Numerical examples 



Example 5.1 First, we consider an artificial example governed by the equation 

^-Au = a(u)\Vu\ A + f, (5.1) 

in the domain f2 = (0, 1) x (0, 1) with cr(u) = 1/(1 +u). The function / is chosen corresponding 
to the exact solution 

u(x, y, t) = 10x(l — x)y(l — y) sech(x + y — t) 2 (5-2) 

which satisfies the homogeneous Dirichlet boundary condition. 

A uniform triangular partition with M + 1 nodes in each direction is used in our computation 
(with h = \/2/M). We solve the system by the proposed method with a linear FE approximation 
up to the time t = 1. To illustrate our error estimates, we take r = h 2 and we present 
numerical results in Tabled! from which we can see that the L 2 errors are proportional to h 2 . 
To demenstrate the unconditional convergence, we take several different spatial meshes with 
M = 16,32,64 for each r = 0.01,0.025,0.05 and we present numerical errors in Table [21 Based 
on our theoretical analysis, in this case, 

\\Uj? -u(;t N )\\ L 2=0(T + h 2 ) \\U?-u(;t N )\\ H l=0(T + h) 

which tend to O(r) as h — > 0. We can observe from Table [2] that for a fixed r, numerical errors 
behave like 0(t) as h/r — > 0, which shows that no time step condition is needed. 

Example 5.2 Secondly, we consider the Burger's equation 

^ + u-Vu-Au = /, (5.3) 

in the unit disk on the plane, with inhomogeneous boundary condition u = g on dVt. The 
functions / and g are given corresponding to the exact solution 

u(x, y, t) = (sech(x + y — t) 2 , cosh(x + y — t) 2 ). (5.4) 

The mesh generated here consists of M boundary points with M = 32, 64, 128, respectively. 
See Figure 1 for the triangulation of the domain. Numerical errors with fixed r and several 
different h are presented in Tables [3] and HI We can see clearly again from Table [3] that the 
numerical errors in L 2 -norm and ff 1 -norm are proportional to 0(h 2 ) and 0(h), respectively, 
when r = 0(h 2 ) and from Table|4]that numerical errors behave like O(r) as h/r — > 0. Thus no 
time-step condition is needed. 

Example 5.3 Finally, we consider the equation 

t^-V ■ (k(u)Vu) = a{u)\Vu\ 4 + f (5.5) 

in $7 = (0, 1) x (0, 1) x (0, 1) with k(u) = 1 + sin 2 u and a(u) = 1/(1 + it). The function / is 
chosen corresponding to the exact solution 

u(x, y, t) = 100x(l — x)y(l — y)z(l — z) sin(x + 2y — z)te~ l . (5-6) 
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Table 1: L 2 -norm errors of the linear Galerkin FEM (Example 5.1). 



T = h' 


h 


\\Uj? -u(;t N )\\ L2 


\\U£ -u{.,t N )\\ m 




1/8 


3.861E-02 


1.657E-01 




1/16 


7.285E-02 


3.211E-02 




1/32 


1.720E-03 


7.678E-03 


convergence rate 


2.08 


2.06 



2: L 2 -norm errors of the linear Galerkin FEM with refined meshes (Example 5.1). 



r = 0.01 


h 


\\U£ -u(;t N )\\ L i 


\\U£ -u(;t N )\\w 




1/16 


9.591E-03 


4.526E-02 




1/32 


5.484E-03 


3.026E-02 




1/64 


4.673E-03 


2.793E-02 


t = 0.025 


h 


\K -u(;t N )\\ L 2 


\\U£ ~u(;t N ) \\ H i 




1/16 


1.569E-02 


8.022E-02 




1/32 


1.167E-02 


6.700E-02 




1/64 


1.079E-02 


6.445E-02 


r = 0.05 


h 


\\U« -u(;t N )\\ L2 


\\Uj? -u(;t N )\\„l 




1/16 


2.486E-02 


1.312E-01 




1/32 


2.079E-02 


1.187E-01 




1/64 


1.984E-02 


1.159E-01 




Figure 1: The FEM meshes with M = 16, M = 32 and M = 64, respectively. 
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Table 3: L 2 -norm errors of the linear Galerkin FEM (Example 5.2). 



1 


M 


||U^-u(-,^)|| L2 


\\U%-u(;t N )\\ H l 




16 


2.138E-02 


1.596E-01 




32 


4.845E-03 


7.534E-02 




64 


1.314E-03 


3.633E-02 


convergence rate 


2.01 


1.06 



Table 4: L 2 -norm errors of the linear Galerkin FEM with refined meshes (Example 5.2). 



r = 0.005 


M 


M -u(;t N )\\ L * 


||U£ -u(;t N )\\ m 




32 


2.112E-02 


1.644E-01 




64 


5.493E-03 


8.883E-02 




128 


4.882E-03 


5.508E-02 


t = 0.010 


M 


-u(;t N )\\ L 2 


||U^-u(.,^)||^ 




32 


2.078E-02 


1.897E-01 




64 


9.551E-03 


1.158E-01 




128 


1.020E-02 


8.863E-02 


t = 0.025 


M 


||U^-u(-,^)|| L2 


||U^-u(.,^)||^ 




32 


2.741E-02 


2.862E-01 




64 


2.520E-02 


2.228E-01 




128 


2.654E-02 


2.051E-01 



which satisfies the homogeneous Dirichlet boundary condition. 

A uniform tetrahedral partition with M+l nodes in each direction is used in our computation 
(with h = y3/M). We solve the system by the proposed method up to the time t = 1. To 
illustrate our error estimates, errors of the numerical solution with r = 8h 2 are presented in 
Table [5j Similalry numerical errors with fixed r and refined h are presented in Table El The 
same observations can be made here. Again, our numerical results show that the scheme is 
unconditionally stabe (convergent). 

6 Conclusion 

We have presented unconditionally optimal error estimates of a class of linearized Galerkin 
FEMs for general nonlinear parabolic equations, which may cover many physical applications. 
The time-step size restriction was always a key issue in previous analysis and practical compu- 
tation. Our theoretical analysis and numerical results show clearly that no time-step condition 
is needed for these linearized Galerkin FEMs. Our approach is based on a priori estimates of 
the numerical solution in the VF 1 ' 00 norm. With these estimates, optimal error estimates can be 
proved unconditionally from classical FEM error analysis. Clearly, our approach is applicable 
to many other time discretization schemes and more general nonlinear equations (systems). 
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Table 5: L 2 -norm rrors of the linear Galerkin FEM (Example 5.3). 



r 


h 


\\U» - U (;t N )\\ L 2 


\\U£? -u(;t N )\\ m 


1/8 


1/8 


2.094E-02 


8.379E-02 


1/32 


1/16 


4.996E-03 


1.983E-02 


1/128 


1/32 


1.220E-03 


4.755E-03 


convergence rate 


2.08 


2.06 



Table 6: L 2 -norm errors of the linear Galerkin FEM with refined meshes (Example 5.3). 



r = 0.025 


h 


\\U£ -u(;t N )\\v 


\\U£ -u{;t N )\\ Hl 




1/8 


1.808E-02 


6.327E-02 




1/16 


4.858E-03 


1.879E-02 




1/32 


1.549E-03 


7.156E-03 


r = 0.05 


h 


\\U» -u(;t N )\\v 


\\U» - U (;t N )\\ m 




1/8 


1.862E-02 


6.707E-02 




1/16 


5.478E-03 


2.357E-02 




1/32 


2.181E-03 


1.206E-02 


T = 0.1 


h 


K-u(;t N )\\v 


\\U» -u(-,t N )\\ m 




1/8 


2.008E-02 


7.762E-02 




1/16 


7.241E-03 


3.759E-02 




1/32 


3.976E-03 


2.668E-02 
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